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ABSTRACT 

The 2001 outburst of WZ Sagittae has shown the most compelling evidence yet for an en- 
hancement of the mass transfer rate from the donor star during a dwarf nova outburst in the 
form of hot-spot brightening. I show that even in this extreme case, the brightening can be 
attributed to tidal heating near the interaction point of an accretion stream with the expanding 
edge of an eccentric accretion disc, with no need at all for an increase in the mass transfer rate. 
Furthermore, I confirm previous suggestions that an increase in mass transfer rate through the 
stream damps any eccentricity in an accretion disc and suppresses the appearance of super- 
humps, in contradiction to observations. Tidal heating is expected to be most significant in 
systems with small mass ratios. It follows that systems like WZ Sagittae - which has a tiny 
mass ratio - are those most likely to show a brightening in the hot-spot region. 

Key words: accretion, accretion discs - binaries: close - novae, cataclysmic variables - stars, 
individual: WZ Sagittae. 



1 INTRODUCTION 

Dwarf novae are cataclysmic variable stars that are observed to un- 
dergo bright outbursts lasting a few days. These outbursts are sepa- 
rated by weeks to months of dim quiescence. At the short orbital pe- 
riod end of the dwarf nova distribution (P or b < 2.2 h), we find the 
SU UMa stars, which also show superoutbursts. These are longer 
outbursts lasting a couple of weeks or more, with the additional 
presence of superhumps, a periodic modulation in the V band light 
curve that rep eats on a time-scale very close to the orbital period. 
IWarnerliT99l gives a comprehensive review of dwarf novae. 

The outbursts of dwarf novae have been explained in terms 
of an instability associated with the gas circulating in the accretion 



disc. Hoshi 1 1979) realised that there is a rapid change in opac- 
ity when hydrogen becomes ionised near 6 500 K. If the temper- 
ature is not higher than this everywhere in the disc, a cyclic be- 
haviour can be established between two thermally stable states: a 
hot, ionised outburst state and a cool, neutral quiescent state. It was 
realised almost immediately that such a thermal change needs to be 
accompanied by a change in the gas viscosity to give the change 
in angular momentum transport that is required to power the ob- 
served accretion rates in dwarf novae. The best candidate for this 
vis cosity remains mag netohydrodynamic turbulence, as suggested 
bv lBalbus & Hawlevl il991l) . Historically, an alternative outburst 
theory was proposed involving an instability in the mass-transfer 
rate from a convectively unstable envelo pe of a heated donor star 
IPaczvnski. Ziolkowski & Zvtkowlll969h . It was the disc instabil- 
ity theory that prevailed due to better compatibility with obser- 
vations, which in nearly all cases have showed no evidence for 
an enhanced mass transfer rate from the companion star. How- 
ever, a re-evaluation of some of the observations iSmakl 19951) and 
evidence from eclipses during the latest outburst of WZ Sagittae 



IPatterson et al 120 02) has led to renewed interest in the possibility 
of enhanced mass transfer. 



In a recent paper, lOsaki & Meved 1200 3l) reviewed the argu- 
ments for and against enhanced mass transfer. They suggested that 
the eclipses observed for WZ Sagittae could be eclipses of the su- 
perhump light source and not the hot spot. The case of WZ Sagit- 
tae is discussed in more detail in Section 2 below, as well as the 
evidence from existing numerical work not discussed by Osaki & 
Meyer. 

lOsakil ll985h originally proposed a model in which super- 
outbursts were caused by an enhanced mass-transfer rate from 
an irradiated seco ndary star. However, numerical experiments by 
White hursJ ll98St) showed that discs in a Roche potential with 
q < 0.3 could spread to a radius at which the material at the edge of 
the disc becomes 3:1 resonant with the tidal field of the secondary 
star. At this point the disc becomes unstable and the gaseous or- 
bits become eccentric. Since superoutbursts are only observed in 
systems with low mass ratios, it seems likely that such a tidal in- 
stability plays a role in their development. Soon after Whitehurst's 
experiments, Osaki formulated the thermal-tidal instability (TTI) 
model lQsakilll989l) . The mass transfer rate from the donor re- 
mains constant in the TTI model, but with every passing normal 
outburst the accretion disc expands until it reaches the 3: 1 resonant 
radius. At this point, Osaki suggested that the rate of angular mo- 
mentum loss becomes more rapid, leading to the accretion of a large 
amount of mass, and a superoutburst. The TTI picture was broadly 
supported by the first two-dimensional (2D) hydrodynamic models 
of superoutbursts iTruss. Murray & WvnifcOOl] ), which confirmed 
that while the outburst is initiated by disc instability in the same 
way as a normal outburst, the tidally-enhanced energy dissipation 
prolongs the outburst, explaining the difference in duration between 
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normal and superoutbursts. A significant eccentricity was generated 
and superhumps did indeed develop on contact with the 3:1 radius. 

In the next section I introduce the evidence for enhanced mass 
transfer in SU UMa stars from observations of hot spot brightness 
and balance this with a discussion of existing theoretical studies of 
mass transfer in these systems. In Section 3 I present the results 
of new, three-dimensional (3D) simulations of tidally unstable ac- 
cretion discs. These show that tidal heating produces a significant 
brightening near the stream-impact region, even with no increase in 
the mass transfer rate through the stream. This means that a simple 
energetic model for hot-spot luminosity is inadequate for a low-q 
system in outburst, as the luminosity in this region is dominated 
by tidal effects. Th e resul ts therefore lend weight to the arguments 
of lOsaki&Meveil fc003l) that the brightening observed near the 
stream-impact region in some SU UMa stars is not due to the hot- 
spot, but is due to the tidal heating of the outer parts of the disc. 



2 ENHANCED MASS TRANSFER: OBSERVATION 
VERSUS SIMULATION 

Until the 2001 superoutburst of WZ Sagittae, little evidence for en- 
hanced mass transfer during dwarf nova outbursts had been found. 
A notable exception is that of VW Hvdrk lVogtl ll983l) found that 
the amplitude of a hump in the light curve increased during normal 
outbursts occurring within 40 days of a subsequent superoutburst. 
This was interpreted as brightening of the hot spot in response to 
an increase in mass transfer rate from the donor star. ISmakl 1 1 9951) 
used similar arguments for an increase in mass transfer rate by a 
factor of two in outbursts of U Geminorum and Z Chamaeleontis, 
employing the relationship 

1 . 2 

ispot = -M2 (VK - V b ) (1) 

where A/2 is the mass transfer rate from the donor star, vk is the 
Keplerian velocity of the gas circulating in the accretion disc and 
Vb is the ballistic velocity of the in-falling gas stream. 

WZ Sagittae is the brightest known dwarf nova in the V-band, 
with V ~ 8 at supermaximum. It is unusual in only showing su- 
peroutbursts, and in its extremely long recurrence time: the 2001 
outburst occurred a mere 23 years after a 33 year cycle of pre- 
vious events in 1913, 1946 and 1978. It has a high-inclination 

(i ~ 75°) and an orbital period of 81.6 minutes iRrzeminski 

1962). The secondary star is of unusually low mass, w ith Mi = 
0.045 ± 0.003 M andg ~ 0.06 iSkidmore et all2002h . The 2001 
outburst was observed in great detail, revealing many layers of vari- 
ability. Superhumps appeared in the light curve 13 days after the 
onset of the outburst, reaching maximum amplitude 24 hours later. 
They persisted for 90 days, well into the decline to quiescence, 
but were preceded by a different m odulati on on the orbital period 
for the first 12 davs. lPatterson et all 12002') found that the depth of 
eclipses in the light curve increased rapidly 6 days after the appear- 
ance of superhumps, and then gradually decreased. This change 
was mirrored by the height of the orbital hump in the light curve. 
The inclination of WZ Sge is not high enough to make eclipses of 
the white dwarf itself observable, so the eclipses in the light curve 
are assumed to be associated with the hot-spot. Patterson et al. con- 
cluded that the hot-spot brightened by a factor 60, which is sugges- 
tive of a commensurate increase in the mass-transfer rate from the 
secondary star. 

I will not discuss the theoretical aspects relating to the re- 
sponse of an irradiated star in terms of mass transfer rate here; 



most r ecently these arguments have been made bv lOsaki & MeveJ 
(2003). Instead, I will discuss briefly the application of analytic 
methods and hyd rodynam i c simu lations to the problem of stream- 
disc interactions. Lubow 1 1994) investigated the response of an 
eccentric accretion disc to an increase in the mass transfer rate 
through a stream. He found that rather than pumping the eccen- 
tricity up further, the addition of low specific angular momentum 
material at a faster rate tended to circularise the gas orbits. Later 
smoothed particle hydrodynamics simulations of dwarf novae with 
mass ratios near the border of tidal stability confirmed that super- 
humps appeared when the mass transfer rate was reduced from 
a rather high value, bu t disappe ared when the rate was restored 
iMurrav. Warner & Wickramasinghell2000l) . These theoretical re- 
sults cannot explain the persistence of the superhumps in WZ Sge 
with an enhanced mass transfer rate. One could circumvent this 
problem if the superoutburst itself was initiated by a burst of mass 
transfer and the disc radius expanded from a value well inside the 
3:1 radius to a radius beyond it. Unfortunately, this argument fails 
for WZ Sge. The observed increase in the eclipse depth occurs six 
days after the appearance of superhumps, so the disc radius has al- 
ready expanded far beyond the 3:1 radius by the time the hot-spot 
region is obse rved to brighten. 

Recently, Foulkes et al. (2004) have used a 2D hydrodynamic 
scheme to model the evolution of spectral line profiles generated by 
an eccentric accretion disc in a binary with q = 0.1. They found sig- 
nificant changes in the profile over the precession cycle of the disc 
and were able to distinguish between the individual components 
from the stream-impact region and the spiral arms. Both compo- 
nents contribute to the superhumps seen in the light curve. 



3 RESULTS 

In this section I present the results from new, three-dimensional 
smoothed particle hydrodynamics (SPH) simulations of the accre- 
tion disc and stream in WZ Sge. The simulations were performed 
on the UK Astrophysical Fluids Facility (UKAFF) supercomputer 
during time awarded in Spring/Summer 2004. 

3.1 Numerical Method 

I use a parallelised version of a 3D SPH code optimised for the 
simulation of accretion discs in compact bi nary sta r s. It is an evo- 
lution of the original code developed by iMurravl ll996h . In the 
SPH method, the fluid equations are solved for each of a set of 
moving points which we call particles. Local physical quantities 
are calculated by means of a weighted sum over all the nearby 
particles within a set distance. This distance, called the smooth- 
ing length, h, is tuned to the local density, making the method ex- 
tremely adaptable in resolving la rge density contras ts. A review of 
the SPH method can be found in Monaghar] 11992}) . In this work, 
h is calculated for each particle such that a minimum of 50 neigh- 
bour particles is maintained regardless of the local density. Here, 
a total particle number of order 400,000 allows a resolved length- 
scale h < 0.017?disc over most of the accretion disc. 

The simulation proceeds by the continuous injection of par- 
ticles from the inner Lagrangian (Li) point. The gas dynamics of 
the particles is calculated in the full 3D Roche potential of a binary 
system, with particles being rejected from the simulation if they 
pass through an inner boundary of the disc at a radius R < 0.05a 
where a is the separation of the two stars, if they intersect the Roche 
lobe of the secondary star or if they achieve the escape velocity at 
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a radius larger than a. The viscosity used in the SP H calculation 
is a refinement of that described by Murrav 1 1996) to mimic an 
Q-viscosity: 

v — ac s H = ac^/Q. (2) 

The lengthscale is chosen to be the scale height H rather 
than the sm oothing length, giving closer compatibility with the 
Shak ura & Sunvaevl Il973l> prescription. 

No viscous calculation is performed for the particles in the 
stream, which move on a purely ballistic orbit. The equation of state 
of all the particles is isothermal, so the rate of energy generation 
in the disc is purely from the viscous dissipation induced by the 
interaction of the particle orbits. Therefore, the effects of radiative 
processes are neglected, making this a simple calculation of the 
total bolometric luminosity due to viscous dissipation at each point 
in the disc. 

3.2 Simulations 

I use parameters representative of WZ Sge: q = 0.07, P or b = 
0.056 days, -M a = 3 x 10~ n M yr" 1 iSmakll 19931 from a 
quiescent hot-spot brightness). The thermal parameters of the par- 
ticles in the disc are chosen as a conservative estimate for a CV disc 
in outburst, that is at the lower end of their probable range: a — 0.1, 
c s = 0.04 ailorb, which assuming a mean molecular mass /i = 0.6 
corresponds to a mid-plane temperature T c = 42, 500 K. The ini- 
tial conditions for t he 3D simulation are generated in the same way 
as in lMurravetalJ i2002). where a 3D disc is replicated from a par- 
ticle distribution derived from an initial 2D simulation. The con- 
ditions are then relaxed to a steady-state before the 3D simulation 
proper. The disc is built up using a constant mass transfer rate of 
3 x 1O _11 M0 yr . The 3D runs start at a point after the outer 
edge of the disc has crossed the 3:1 radius and superhumps have 
developed. Three runs are then performed with initial conditions 
that are identical except for the rate of mass injection. Run A con- 
tinues with the same mass transfer rate. Run B proceeds with a rate 
five times lower, run C with a rate five times higher. The number of 
particles at the start of each run was N = 428685. 

Fig0 shows the evolution of the viscous energy dissipation 
rate and the (k,l) = (1,0) eccent ric mode, calculated using the 
method described in Section 3.3 of iTruss. Murrav & Wvnnl 1200 ll 
hereafter TMW). The superhumps are clearly visible in the run with 
the unchanged mass transfer rate (A) and in the run with the re- 
duced rate (B). However, they disappear when the rate is increased 
in run C, with the eccentricity falling throughout the simulation as 
the disc edge is driven away from the 3: 1 radius. The difference in 
length of the simulations is due to the finite supercomputer time. 

FiglH shows a map of the energy dissipation rate in the disc 
and a corresponding Doppler velocity map at time t = 8 days 
for run A. The hot-spot region and tidally-induced spiral struc- 
ture are clearly evident, as are velocity components consistent with 
gas on ballistic and Keplerian trajectories. There is also a non- 
Ke plerian component i n the velocity map, similar to that noticed 
by Foulkes et al. (2004). The radial luminosity profile of the disc 
is presented in Fig|3|for the two sectors shown in the dissipation 
map. This has been calculated by dividing the disc up into annu- 
lar bins of width 0.01a and averaging the luminosity (dissipation 
rate per unit area multiplied by the area of each bin) in the az- 
imuthal direction. Also plotted is the expected total bolometric hot- 
spot luminosity calculated from equationQfor three different mass 
transfer rates. Even with this small emitting area per radial bin, the 
luminosity in the hot spot region is well above that expected for 
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Figure 1. Viscous energy dissipation rate per unit mass and eccentric mode 
strength during the three simulations. A: — M% = 3 X 1O~ 11 M0 yr — 1 , 
B: -M 2 = 6 X lCr 12 M Q yr" 1 , C: -M 2 = 1.5 X 10- 10 M Q yr" 1 . 



—Mi = 3 x 10 _11 Mq yr -1 . If one assumes that the radial extent 
of the hot-spot covers the region of tidally-enhanced dissipation be- 
tween R = 0.5a and R = 0.6a, the total 'hot-spot' luminosity is 
nearly an order of magnitude higher still. 

A feature that becomes apparent from the 3D simulations is 
the magnitude of the eccentricity, which is far smaller than has 
been found in previous 2D work. The mode strength in FigQis 
of the order a few times 10 -3 , in comparison with ~ 0.5 for 2D 
superoutburst simulations of Z Cha (TMW). As a result, the pattern 
of dissipation shown in Fig.|2|does not change appreciably over the 
course of the simulation, and the dissipation remains enhanced in 
the region of the hotspot throughout. However, the FigQshows that 
the eccentricity of the 3D discs in runs A and B are increasing at 
the end of the simulations. It remains unclear whether this trend 
would continue if the simulations were allowed to run indefinitely, 
although the growth rate of the eccentric mode in 3D does not sug- 
gest that the eccentricity will reach the levels achieved in the 2D 
case over a sensible time-scale for a dwarf nova outburst. 

Despite the low eccentricity, a clear superhump signal is still 
seen and the tidal interaction still has a significant influence on the 
luminosity generated in the hot-spot region. Of course, if the ec- 
centricity were higher the interaction would be even stronger. We 
would also expect to see more modula tion in the diss ipation map 
over a precession cycle, as discussed by Foulke s"et alj 120041) . 

As a final check, a short simulation was performed for a disc 
not subject to the tidal instability, with an outer radius well in- 
side the 3:1 resonance. The results are shown in Figs. [4] and [5] 
In this case, the structure of the hot spot is exactly as expected. 
There are clearly defined components of emission at the ballistic 
and Keplerian velocities and very little else, only a very weak spiral 
structure. The integrated luminosity of the hot-spot is now consis- 
tent with equation Q for the mass transfer rate that has been mea- 
sured in quiescence, and that is used in the simulation: —Mi = 
3 x lO-^Moyr -1 . 



4 M. R. Truss 




Figure 2. Map of energy dissipation rate in the simulation with — M2 = 3 X 10 — 11 Mq yr" 1 and corresponding Doppler velocity map. The straight lines 
on the left define the sectors that are considered in generating the radial profiles in Fig.|3] The Roche lobe of the secondary star and velocities of the ballistic 
stream and its Keplerian equivalent are plotted on the right for q = 0.07. 



4 DISCUSSION 

I have shown that tidal heating accounts for a significant fraction of 
the energy dissipated near the stream-impact region in an eccentric 
accretion disc. The tidal heating component dominates that due to 
the impact of the ballistic stream, producing an apparent brighten- 
ing near the hot-spot. No increase in the mass transfer rate from the 
secondary star is required to produce this brightening. Therefore, 
even for a disc with low eccentricity the simple model for the lumi- 
nosity of the hot-spot given in equationQis inadequate to estimate 
the mass transfer rate from the donor. Of course, this equation is 
not incorrect: but for a small correction due to the slightly non- Ke- 
plerian velocity of the orbiting gas in the disc, the luminosity of the 
impact alone will be given correctly. It is j ust that this lum inosity 
is swamped by the tidal heating component. Osaki & Mever 1 2003) 
suggested that it is the eclipse of the superhump light source (dom- 
inated by the tides) and not the hot-spot that has been observed in 
WZ Sge. This work supports this view, with the caveat that it could 
well be that the eclipse is that of the hot-spot, the luminosity in this 
region is simply being dominated by the tides. For a disc with zero 
eccentricity, such as a quiescent disc with an outer radius much 
smaller than the 3:1 radius, the expression in equation Qremains a 
good diagnostic of the mass transfer rate. 

This brightening of the hot spot region will be observed in 
systems with mass ratio q — M2 jM\ < 0.3, where tidal instability 
becomes important. It may seem intuitive to take this assertion a 
step further; that the lower the mass ratio, the stronger the tidal 
heating and the brighter the hot-spot. However, it should be stressed 
that to date no systematic studies of the variation in strength of the 
tidal heating component with mass ratio have been carried out. 

In terms of the sequence of events during an outburst of WZ 
Sge, the picture that has emerged from this work would be as fol- 
lows. I assume that the outburst is initiated in the normal way by 
the disc instability, with the disc spreading past the 3:1 radius some 
time later. The growth rate of the eccentricity must be very fast 
as the superhumps were observed to reach maxi mum amp litude 
within 24 hours of appearing in the light curve IPatterson et alj 
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Figure 3. Radial profile of luminosity per radial bin of width 0.01a for the 
two sectors in Fig. [2] The solid line has been smoothed over the sector that 
includes the hot-spot, the dashed line is the equivalent for the diametrically 
opposite sector. The numbered curves are the expected total bolometric hot- 
spot luminosity calculated from equation0for — M2 = 3 X 10" 11 (1), 
10" 10 (2) and 3 X 10- 10 M Q yr" 1 (3). 

2002). The brightening effect due to the tides will be first observed 
up to one precession cycle later when the apastron of the disc passes 
L\. This is consistent with the delay of a few days that was ob- 
served between the appearance of the superhumps and the bright- 
ening event. As the disc empties and the outer radius finally shrinks 
inside the 3:1 radius, the effect will diminish over the course of 
the outburst until the brightness of the hot-spot becomes consis- 
tent with the value given by equation Q for the mass transfer rate 
measured in quiescence. 

I have confirmed with a 3D model the analytic 
predic tion of iLubowj tl994|) an d the 2D prediction of 
Murrav. Warner & Wickramasinghe (2000) that an increase 
in mass transfer rate from the donor drives the disc away from the 
3:1 resonance and suppresses the appearance of superhumps. This 
is in contradiction to what has been observed in WZ Sge. There are 
interesting differences between the 3D results presented here and 
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Figure 4. As Figl21 but for a smaller disc with -Rdiac < ^?3:i- The ballistic and Keplerian components are much more clearly defined in the Doppler map and 
there are only very weak spiral features. The straight lines on the left hand side define the sector of the disc used to calculate the luminosity profile in Fig|5l 
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Figure 5. As Fig|3] but for a smaller disc. The integrated hot spot luminosity 
is consistent with the curve calculated from equation with — M.% = 3 X 

lO^Mgyr" 1 . 



previous 2D work. In particular, while eccentricity and superhumps 
develop in the same way in 2D and 3D, it is much more difficult to 
generate high eccentricities in 3D for reasonable values of sound 
speed and viscosity parameter. Indeed, the eccentricity of e ^ 0.3 
estimated for WZ Sge by Patt erson et alJ l2002h is not found in 
these 3D simulations. However, the tidal heating effect is still 
significant even with a small eccentricity. If future observations 
can reveal with accuracy the exact shape of the disc and e is found 
to be as high as 0.3, we will need to consider how this can be 
achieved. A possible solution may lie in a magnetic field anchored 
on the rotatin g white dwarf, which could have a pr opelloring effect 
on th e gas iLivio & Pringld 1 1994k lLasota. Kuulkers & Charlesl 
1999) that keeps the outer edge of the disc past the 3:1 radius 
and allows a large eccentricity to develop. I will investigate this 
numerically in a forthcoming paper. 

As an aside to these issues, it is interesting to note that an 
eccentric disc has no well-defined single outer radius. Indeed, if we 
assume the disc is an ellipse with an eccentricity e, the ratio of the 
semi-major to semi-minor axes is 



So, for e ~ 0.3, there is nearly a factor of 2 variation in apparent 
radius over a single precession cycle. This should always be borne 
in mind when inferring accretion disc radii from observations in 
tidally unstable systems, especially given that the precession pe- 
riod can be several days in comparison to the orbital or superhump 
periods, which are usually shorter than a couple of hours. 

In conclusion, it is apparent from this work that the observa- 
tions of hot spot brightnesses during dwarf nova outbursts can be 
attributed to the interaction of a constant- M accretion stream with 
the edge of an eccentric accretion disc. Future work on the theo- 
retical side should take into account the effects of radiative trans- 
fer in the spot region. This will be a challenge, and has not been 
considered in this work, however it is expected to have an impact 
on the predicted brightnesses. In particular, the amount of stream 
overflow above and below the rim of the disc is sensitive to the ra- 
diative cooling rate in the hot-spo t region lArmitage & Livid 1998t 
iKunze. Speith & Hessmanj2 000). This will have observable conse- 
quences on the shape of the measured eclipses and the appearance 
of the orbital hump in the light curve. From an observational point 
of view, it is desirable to discover the exact dimensions and geome- 
try of a precessing disc during outburst. This will be a robust test of 
the TTI model for superoutbursts, which relies on the disc edge en- 
countering the 3:1 resonance and the mass transfer rate remaining 
fixed. 
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